Disordered ground states in a quantum frustrated spin chain with side chains 
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We study a frustrated mixed spin chain with side chains, where the spin species and the exchange 
interactions are spatially varied. A nonlinear a model method is formulated for this model, and a 
phase diagram with two disordered spin-gap phases is obtained for typical cases. Among them we 
examine the case with a main chain consisting of an alternating array of spin-1 and spin-| sites 
and side chains each of a single spin-| site in great detail. Based on numerical, perturbational, and 
variational approaches, we propose a singlet cluster solid picture for each phase, where the ground 
state is expressed as a tensor product of local singlet states. 
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I. INTRODUCTION 

Quantum one-dimensional (ID) spin systems have 
been studied in various aspects, especially with interest 
on their strong quantum fluctuations due to the low di- 
mensionality. There appear a variety of quantum disor- 
dered ground states where the continuous spin rotation 
symmetry is not broken and the lowest spin excitation 
has a finite gap (spin-gap). These quantum disordered 
states have no analogues in classical spin systems. Typi- 
cal examples are a Haldane state in a spin-1 chain [1], a 
dimer state in a spin-^ chain with bond alternation [2], 
and a spin-gap state in a spin-^ ladder [3]. 

In extensive research for various ID spin systems, spin 
chains with side chains have not attracted enough atten- 
tion in spite of its potentially rich physics. Since a side 
chain is of finite length, it may enhance quantum fluc- 
tuation in the system. Actually, the ID Kondo necklace 
model, which has been extensively studied as a simplified 
version of ID Kondo lattice model [4, 5], can be regarded 
as a spin chain with side chains. In this model, the main 
chain is a spin-| chain and each side chain consists of a 
single spin with magnitude ^. The ground state of this 
model is known to be in the Kondo singlet phase with 
spin-gap [5], while the spin-^ chain without side chains 
is critical. This means that the quasi-long range order in 
the main chain is destroyed by the quantum fluctuation 
in the side chains. Also, if the side chains bring geomet- 
rical frustration into the system, quantum fluctuation is 
expected to be further enhanced. Thus, it is an inter- 
esting subject how quantum fluctuation manifests itself 
and what kind of ground state appears in various types 
of spin chains with side chains. 

In this paper, we investigate the natures of quantum 
disordered ground states of one of the simplest models 
with frustrated side chains. The main chain of the model 
consists of two species of spins in alternating order, and 
each side chain consists of a single spin which are alter- 
nately attached to the main chain, as will be shown in 



Fig. 1 of the next section. This model incorporates the 
effects of mixture of different spins, bond alternation, 
and frustration in spite of its simplicity. In particular, 
the frustration comes from triangles, each consisting of 
antiferromagnetically interacting three spins. Although 
quantum spin systems with similar geometry have been 
investigated by several authors [6], our model is physi- 
cally different from them. If the side chains of our model 
are removed, the main chain is in a ferrimagnetic ground 
state. The side chain spins introduce frustration to the 
system, and destroy the ferrimagnetic order leading to 
quantum disordered states. 

The disordered ground states of the present model can- 
not be understood in the conventional valence bond solid 
(VBS) picture [7], which successfully explains the dis- 
ordered ground states of many spin models with local 
frustration. Instead, we will explain the present ground 
states in the concept of the singlet cluster solid (SCS) 
picture. A SCS state is a direct product of local singlet 
states, or singlet clusters. Each singlet cluster consists 
of more than two singlet dimers, and the dimers are res- 
onating locally within the cluster. The SCS states are 
realized as a result of the interplay of quantum fluctua- 
tion and local frustration, as will be explained in detail. 

It is desirable that the SCS state manifested in this pa- 




FIG. 1: A quantum spin chain with side chains; two unit cells 
are presented. Si, S2, and T are spins whose magnitudes 
satisfy Eq. (3). The case that Si — 1^ 82 — \^ and T = ^ is 
studied particularly in detail. 
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per is experimentally inspected in materials. However, a 
material precisely described by the present model is not 
found so far within our knowledge. Despite the lack of 
materials, it is worth clarifying the concept of the SCS 
states and verifying the existence of them in a concrete 
model. Further, considering the rich variety of magnetic 
materials synthesized by the modern chemical technol- 
ogy [8], desired materials are expected to be synthesized, 
since they are not necessarily complex in structure. 

The details of the model are explained in the next 
section. Various approaches are employed to clarify the 
ground state phases of this model: In section III, a non- 
linear <j model (NLSM) method is proposed to grasp 
qualitative feature of the phase diagram. Since a sim- 
ilar NLSM method has been developed for mixed spin 
chains without side chains so far, we extend it to the 
present side-chain case. 

We also employ other approaches limiting ourselves to 
the simplest case of spin magnitudes 1 and ^. In section 
IV, the numerical diagonalization for finite systems is 
carried out to obtain the quantitatively reliable phase di- 
agram, which confirms the qualitative correctness of the 
NLSM method. Some limiting cases are exactly treated 
to draw physical picture for each phase in section V. Af- 
ter these considerations, we arrive at the SCS pictures 
to explain the ground states in section VI. The SCS 
pictures are also supported by variational calculations in 
section VII. The last section is devoted to summary and 
discussion. 



II. MODEL HAMILTONIAN AND ITS 
CLASSICAL VERSION 

We study an isotropic quantum spin chain with al- 
ternately arrayed side chains as is illustrated in Fig. 1. 
In the pth unit cell, Si{p) and S2{p) are spin operators 
on the main chain, and T{p) is a spin operator on the 
side chain. The quantum numbers of spin magnitudes 
of these spin operators are expressed as ^i, 5*2, and T, 
respectively. Exchange parameters are represented as Ji , 
J2, Ki and i^2, and assumed to be all positive. Then the 
Hamiltonian is written as 

N 

H = J2{Ji Mp) • Mp) + J2 S2{p) • Si(p + 1) 

+ K, Si(p) • T{p) + K2 S2{p) • T{p)} . (1) 

The spacing between nearest spins is a, and the length 
of a unit ceU is 2a. The Hamiltonian is characterized by 
three independent dimensionless parameters: 



J 



k 



TKi 



S2K2 



(2) 



Here j measures the strength of the bond alternation in 
the main chain, k measures the strength of interaction 
between a main-chain spin and a side-chain spin, and r 
measures the strength of frustration [9]. 



In the present paper, we assume the following restric- 
tion on spin magnitudes: 



Si-S2-T = 0. 



(3) 



This is the condition that the corresponding classical spin 
chain can have a ground state with no total magnetiza- 
tion (i. e. no ferrimagnetism) , when K2 is not large. The 
restriction (3) serves to simplify the Berry phase term in 
the continuum limit. 

Expectation values of the spin operators for a spin co- 
herent state are written as 



(Si(p)) = 5iMi(p), 
(S2(p)) = -52M2(p), 
(T(p)) = -TM^(p), 



(4) 



where Mi(p), M2(p), and M.±{p) are unit vectors. Re- 
placing the spin operators in Eq. (1) by them, we have 
the classical version of the Hamiltonian: 



N 



= [^2(P) - Mi(p + 1)]2 



+ 1 



+ 



kr 
1 -r 



[Mi(p)-M2(p)]^ 



^ [Mi(p) - rM2(p) - (1 - r)m^{p)f], 



(5) 



where Ji = J1S1S2 and a constant term is omitted from 
He- The classical antiferromagnetic configuration, 



Mi{p)=M2{p') = M^{p") 



(6) 



for all j9, y, and p'\ is the ground-state solution if the 
pre-factors of the squares in Eq. (5) are all positive. This 
gives the condition for the classical stability of the anti- 
ferromagnetism as 



0<r<l, 0<A;<--1. 

r 



(7) 



In the following arguments, we will concentrate on this 
region. 



III. NONLINEAR a MODEL FOR THE SPIN 
CHAIN 

Using the spin coherent representation, the partition 
function of Hamiltonian (1) is written in a path-integral 
form as 



Z = lY[D[Mj]S{M^j-l)e- 



(8) 
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with j = 1,2, and _L. The action A is written as 

A = -i^B + Ah, 

= Y,{Siw[M,{p)] - S2w[M2{p)] - Tw[M^{p)]}, 



Ah 



-f 

Jo 



drHc, 



(9) 



where f3 is the inverse of temperature, and w[M.i{p)] is 
the sohd angle which Mi(j9) forms in period p. The term 
—iAb in action A is the Berry phase term. 

We introduce slow variable m{p) for each unit cell, and 
fluctuation variables Li(p), L2(p), and L^(p) for each 
spin in a unit cell. Then the original variables are trans- 
formed as follows: 

Mi(p) = m(p) + aLi(p), 

M2(p) = -m{p) + -m(2^+ 1) + aL2(p), 



M_l(p) =m{p)^aL^{p). 



(10) 



This transformation is found by observing and extend- 
ing the transformation for a simple chain without side 
chain [11, 12]. Since the left hand sides in (10) are unit 
vectors, we have the following constraints for the new 
variables: 

= 1, mLi=mL2 = m- = 0. (11) 

The fluctuation variables depends on one another, and 
one of them, e. g. L_l can be set equal to 0. Hence 
the number of independent variables are conserved in the 
transformation. Deflning new fluctuation variables. 



R = L2 — Li, 

Q = (l-r)Li-Li+rL2, 
and taking the continuum limit, we have 

A = jdTdxi^^S2^^ ■ (too. X 



(12) 



+ -Ji(c+-A:r^) 



dx 



a - 



Ji[kW + 2f • Q) + c+(R^ + 2g • R)] j (13) 

with c± = 1 =b j — A:r/(1 — r). Here vectors f and g are 
given as 



_ dm . Si — S2 dm 
I = r— h I— m X — — , 

ox aJik{l - r) or 

C- dm . rSi — S2 dm 
g = I m X — — . 

c+ dx c+aJi(l - r) dr 



(14) 



Integrating the partition function with respect to R 
and Q, we have the following NLSM action: 



leff 



dx{ i — m • 

47r 



dm dm 

dr dx 



c+(5i - ^2)' + k{rSi - 82)^ [dm 



+ - Jic+ ( 1 



4aJi(l - rykcj^ 
a - ( ^ (?- \ ( dm\^ 



dr 



(15) 



The first term is the topological term and Q is the topo- 
logical angle given by 



47rj(52-r-5i) 
(1 - r)(l + j) - kr' 



(16) 



In the absence of frustration (r = 0), the topological 
angle reduces to 6> = 47r5'2J2/(Ji + -^2)- This is the same 
expression as that for a simple spin chain of magnitude ^2 
with bond alternation [10, 11]. This can be interpreted 
as follows: a spin Si is combined with the adjacent spin 
T on the side chain, and a two-spin cluster with total 
spin magnitude S\ —T (= 5*2) is formed. However, the 
coeflicients of {dm/dr)'^ and {dm/dx)'^ in Eq. (15) do 
not reduce to those for the simple spin chain, even if k is 
very large. This reflects that the quantum fluctuation of 
spins on side chains still survive for large k. 

The topological term of Eq. (15) determines whether 
or not the system has a spin-gap in the same manner as 
Haldane's argument [1]. That is, the system does not 
have a spin-gap if and only if the topological angle is 
just TT (mod 27r). This condition is written as 



2j(^2-r^i) 
(1 - r)(l + j) - kr 



(17) 



where h is any half odd integer. This gapless condition 
determines phase boundaries between gapful disordered 
phases in the parameter space. We notice that, for each 
value of /i, boundaries for all values of r pass through the 
common point 



2{Si - S2) h 



252-/1 '252-/1 



(18) 



In the case of Si = 1^ S2 = ^ and T = |, only a 
permitted value of h in Eq. (17) is ^ for /c > and 
j > 0. Then the phase boundaries for several values 
of r are solid lines in Fig. 2. Owing to the deflnitions 
of the parameters, they are straight lines in the present 
approximation. The regions of the both sides of each 
boundary are gapful disordered phases. We call them 
Gap I phase and Gap II phase as noted in the flgure. 
For r = the phase boundary is horizontal, since the 
topological angle is independent of k as mentioned below 
Eq. (16). This is understandable by considering Si as a 
composite of two ^ spins. In fact, when r = 0, T and 
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IV. NUMERICAL DIAGONALIZATION 




FIG. 2: Phase boundaries for several values of r in the k-j 
plane by the NLSM method for 5'i = 1, 5*2 = |, and T = |. 
For each value of r (= K2/2K1) the region ofO</c<^ — 1 
is meaningful in the NLSM method. 



one of the | spins of Si necessarily forms a valence bond 
irrespective of the value of k in the ground state. As 
r increases from 0, the slope of the boundary becomes 
negatively large. The reason will be argued later. 



We also show phase diagrams in other cases in Fig. 3; 
the phase boundaries in (a) are for = |, ^2 = 1, and 
T = |, and those in (b) are for S*! = |, 6*2 = ^, and 
T = 1. Although the equation (17) determining phase 
boundaries are quite general, we mainly examine the case 
of S^i = 1, ^2 = |, and T = |. This case is expected 
to include the essence of the present type of spin chains 
with side chains. 
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Hamiltonian (1) can be numerically diagonalized for 
small size systems. The numerical calculation is effective 
not only to analyze the system itself in detail, but also to 
know the preciseness of the NLSM method by comparing 
the results. We performed numerical diagonalization in 
the case of 5'i = 1, 5*2 = ^, and T = | to obtain the 
phase diagram for the ground state. 

The phase transition points are determined as follows. 
The phase transitions between different spin gap phases 
are expected to be the Gaussian transition. Hence we 
employ the method of twist boundary condition proposed 
by Kitazawa [13] and Kitazawa and Nomura [14] to de- 
termine the phase boundary. As will be examined later, 
the ground state phases are described by different SCS 
configurations. Under the twisted boundary condition, 
the different singlet solid configurations have different 
time reversal parities depending on the even-odd parity 
of the number of valence bonds across the twisted bound- 
ary. Hence the energy levels of the ground state and the 
first excited state cross at the phase boundary without 
level repulsion. This ensures the precise evaluation of the 
phase boundary. The size extrapolation is based on the 
following formula for the finite size correction [13, 14]: 



ic(00)=ie(iV) + |, + ^, 



(19) 



where jc{N) is the finite size critical value of quantity j, 
and ci and C2 are fitting parameters. We have carried 
out the extrapolation using numerical results for total 
spin number SAT = 12, 18, and 24. 



0.0 0.5 1.0 1.5 2.0 
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FIG. 3: Phase boundaries in the k-j plane by the NLSM 
method in the cases of (a) 5*1 = |, & = 1, and T = |, and 
(b) in the case of 5*1 = |, 6*2 = |, and T = 1. For each value 
of r the region ofO</c<^ — lis meaningful. 



FIG. 4: Phase boundaries by numerical diagonalization for 
Si — 1^ S2 — and T — ^. Each point of a boundary is 
determined by extrapolation for the total spin number 3A/' = 
12, 18, and 24. Inset: Phase boundaries for small k. 



5 



Resultant phase boundaries for several values of r are 
plotted in Fig. 4. Comparing Figs. 2 and 4, we find 
that the NLSM method gives qualitatively correct phase 
boundaries. In particular, for small r, or weak frustra- 
tion, the NLSM method provides a quantitatively fair 
approximation. With the increase of r, the Gap I phase 
extends to the region j > 1 for small /c, and it is sup- 
pressed to the region j < 1 for large k. Although this 
feature qualitatively coincides with that of the NLSM re- 
sults, quantitative coincidence becomes worse with the 
increase of r. This is natural because the present NLSM 
method starts from a classical antiferromagnetic solution 
in the absence of frustration. 

The possibility of the first order transition between dif- 
ferent spin-gap phases has been pointed out in the frus- 
trated ladder by Hakobyan and coworkers [15] in the ap- 
propriate parameter regime. Considering the presence of 
frustration, this type of transition cannot be ruled out 
in the present model. However, we did not find the nu- 
merical evidence for the first order transition within the 
parameter regime discussed in this paper. 



LIMITING CASES 



(a) 




Si J-\ S2 J2 



(b) 




Si Ji S2 J2 



FIG. 5: (a) 3-spin units for j 0, and (b) those for j ^ 00 
and r ^ in the case of 5*1 = 1, 6*2 = ^, and T = |. Three 
spins (circles) in a shadowed region form a 3-spin unit. 



This is a singlet eigenstate of Hs belonging to the eigen- 
value 



Eq = —Ji 



(22) 



To further confirm the numerical phase diagram for 
Si = 1, S2 = and T = ^, we consider the effective 
theory in the limiting cases of j ^ and j ^ 00. 



A. Strong Ji limit {j 0) 



(ii) For = 1, we have = 1, 0, or —1. In this 
three-dimensional subspace, it is sufficient to inspect the 
case of = 1 owing to the spherical symmetry of the 
Hamiltonian Hs. We choose the orthonormal basis of the 
subspace 5* = -S^ = 1 as 

||i,i,i)) = ^(|^iT)+I^Ti)-V2|0TT)), 



In the limit of j ^ 0, the system can be regarded as a 
one-dimensional array of weakly coupled 3-spin units as 
shown in Fig. 5(a). One of the 3-spin units is described 
by the Hamiltonian 



Hs = JiSi • S2 + i^iSi • T + K2S2 ■ T, 



(20) 



where we have dropped the common index p representing 
the pth unit cell for simplicity. The Neel basis is repre- 
sented as I 5^ , , ), where Si takes or J|, and 
and T take t or |. By introducing the composed spin 
S = Si + S2 + T, we have another set of basis vectors 
II 6*, 5*^, a)), where S and S^ are quantum numbers of the 
magnitude and the z-component of S respectively, and a 
discriminates multiple states with the same S and S^, if 
necessary. We seek the ground state of the 3-spin unit 
with = or 1, since all the exchange interactions are 
antiferromagnetic . 

(i) For S' = 0, we have S^ = 0. Then the one- 
dimensional subspace consists of a single state. 



io,o)) = -^(ioiT)+ion; 



V2\ifii)-V2\m))- 



(21) 



iii,i,2)) = — (itiT)- mi))- 



(23) 



Operating Hs on these bases, we have an eigenvalue equa- 
tion. Then the lowest eigenvalue Ei in this subspace is 
determined as the smaller solution of the characteristic 
equation: 



0. (24) 



Introducing a normalized energy difference as e 
{El - Eo)/Ji, Eq. (24) with Eq. (22) reduces to 



2e' 



[3(1 + A:) - 4rk]e + 2A:[2 - (1 + A:) r] = 0. (25) 



If the smaller solution for e is negative, the ground state 
is a triplet (.9 = 1) state; otherwise it is a singlet {S = 0) 
state. Using (25), the condition for the triplet ground 
state becomes 



k ^ kn 



'--1 

r 



(5 = 1). 



(26) 



We notice that kc > 1 in the region of < r < 1, which 
we have concentrated on in this paper. In the triplet 
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ground state, S2 tends to orient the opposite direction to 
Si for Ji > i^i, and T does for Ji < Ki. The composed 
spin S always orients to the same direction as Si. 

The composed spin S(p) at the pth unit cell interacts 
with adjacent S(p + 1) by an effective exchange interac- 
tion. We denote the effective exchange parameter by Jeff. 
Since the interaction between S2(p) and Si(p + 1) is an- 
tiferromagnetic (J2 > 0), the correlation between Si(p) 
and Si(p + 1) is antiferromagnetic for i^i > Ji (A: > 1) 
and ferromagnetic for Ki < Ji {k < 1). Therefore Jeff 
has the same sign as Ki — Ji = Ji(/c — 1), considering 
the signs of (S(p) • S{p + 1)) and (Si(p) • Si(p + 1)) are 
the same. 

For k > /cc, we have Jeff > 0, since /cc > 1 for < 
r < 1. Hence the original spin chain is equivalent to a 
spin-1 antiferromagnetic Heisenberg chain consisting of 
effective spins, S(p)'s. The ground state of an uniform 
spin-1 chain is the Haldane state [1], which gives a spin- 
gap for excitation. In the Haldane state, there is strong 
correlation on each adjacent spin pair, as known from 
the VBS picture for effective spins, S(p)'s [7]. In terms 
of the original spins, there is strong correlation between 
adjacent 3-spin units. For /c < /Cc, on the other hand, the 
ground state of each 3-spin unit is already a closed singlet 
state. Then the ground state of the total spin chain is 
approximately an array of such closed local singlets, and 
there is almost no correlation between adjacent 3-spin 
units. Thus there is a Gaussian transition between the 
two characteristic ground states with spin-gap dX k = kc- 
The value of kc in this argument for j ^ agrees with 
the critical value by the numerical diagonalization, as is 
seen on the j = line of the phase diagram (Fig. 4). 



B. Strong J2 limit {j 00) 

In the limit of j 00, spins Si(p+ 1) and S2(p) form 
an effective spin S(p) = Si(p + 1) + ^2{p) with mag- 
nitude |, and other interactions can be treated as per- 
turbations. Then the effective Hamiltonian for S(p) and 
T{p) is 

ffeff = ^|-gJlS(p) -8(^+1) 

+ ^ifi S(p) • T{p + 1) - ^-K2 S{p) ■ T(p) I . (27) 

The ground state of this chain is still nontrivial. However, 
K2 plays a secondary role in the weakly frustrated region, 
so that each T antiferromagnetically interacts with the 
ferromagnetic chain consisting of S's. Therefore, it is 
plausible that the ground state is always nonmagnetic 
for small r. Numerical studies of the effective model (27) 
suggest no phase transition for < r < 1/3 where no 
phase transition is predicted by the NLSM method for 
large j. 
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FIG. 6: Four regions in a typical phase diagram for 6*1 = 1, 
& = ^, and T = ^; the phase diagram for r = 0.2 by the 
NLSM method is shown. The solid line is the phase boundary 
between Gap I phase and Gap II phase. In each phase, the 
dashed line oi j — 1 means a crossover between regions with 
different features. Gap I phase consists of regions I-A and I- 
B, while Gap II phase consists of regions II- A and II-B. They 
are explained by SCS pictures. 

In terms of the original Hamiltonian (1), the ground 
state of j ^ 00 and r ^ is a direct product of local 
singlet states of 3-spin units. A 3-spin unit consists of 
S2(p), Si(p+1), andT(j9+l), as shown in Fig. 5(b). The 
ground state of finite k and r (j ^ 00) is adiabatically 
connected to the limit without phase transition as long as 
r is small. The full phase diagram of the effective model 
(27) is investigated in a separate paper [18]. 

VI. SINGLET CLUSTER SOLID PICTURE 

In this section, we propose the SCS picture to ex- 
plain any ground state in the phase diagram for Si = 1^ 
S2 = \^ and T = |. This is a generalization of the VBS 
picture, and is based on expressing Si with magnitude 1 
as 

Si = S^'^+sf\ (28) 

where S^^^ and S^^ are spins with magnitude ^ [16]. 

For convenience of explanation, we divide each phase 
into two regions by the line of j = 1 as schematically 
shown in Fig. 6: Gap I phase is divided into regions I-A 
and I-B, and Gap II phase is divided into regions I-A and 
I-B. 



A. VBS picture and its insufficiency 

The ground states of the limiting cases in the preceding 
section are explained by VBS pictures. The VBS picture 
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FIG. 7: VBS pictures in the limiting cases for Si = 1, S2 = 
and T = |: (a) VBS I represents the ground state for j ^ 
and k < kc, and (b) VBS II does for j ^ 00. Small circles 
are |-spins, and bold gray lines are valence bonds. Spin Si is 

expressed by two |-spins as Si = S^^^ +S^^\ Loops including 
two valence bonds are for correspondence to SCS pictures (see 
text). 



for j ^ and A: < kc is VBS I illustrated in Fig. 7(a) [17]. 
The valence bonds on Ji -interactions, which we hereafter 
abbreviate as the Ji -valence-bonds, mainly contribute to 
the energy gain of the ground state of VBS 1. On the other 
hand, the VBS picture for j ^ oc is VBS II illustrated in 
Fig. 7(b). The valence bonds on J2 -interactions, or the 
J2 -valence- bonds, mainly contribute to the energy gain 
of the ground state of VBS 11. 

The VBS pictures are not adequate for regimes away 
from the above limiting cases, although they are expected 
to be qualitatively valid for small- j and large- j regimes. 
As seen in Fig. 6, the energetic advantage of Ji-valence- 
bonds of VBS I on line j = {k < kc) remains within 
region I- A because Ji > J2. However the advantage is 
lost in region I-B because J2 > Ji- Similarly, the ener- 
getic advantage of J2- valence- bonds of VBS II in the limit 
of j ^ 00 remains within region II- A because J2 > Ji. 
However the advantage is lost in region II-B because 
Ji > J2. Since the line of j = 1 is not a phase bound- 
ary, we need a new picture to explain the whole Gap I 
(II) phase which reduces to VBS I (II) in the limit. The 
picture will be a SCS picture. 



B. Concept of SCS picture 



FIG. 8: (a) SCS I, the SCS picture for Gap I phase {Si = 1, 
5*2 = |, and T = |). SCS I is a tensor product form of 
singlet clusters, (b) A singlet cluster \ipi{c)) (the left hand 
side) in SCS 1. It is represented as a linear combination of two 
valence bond states, | ai ) and | pi ) (the right hand side), c is 
the coefficient of the linear combination. 



Here we have considered that any spin with magnitude 
more than 1 is resolved into a set of spins with magni- 
tude ^. Then a singlet cluster is a singlet state of more 
than two spins with magnitude |. 

A singlet cluster in a SCS state is represented as a 
superposition of products of valence bonds. Hence the 
valence bonds are resonating within the singlet cluster. A 
VBS is a special case of the SCS, where a singlet cluster 
is a single valence bond and no resonation occurs. A 
resonating valence bond (RVB) state is another special 
case, where the whole system is the singlet cluster and 
all valence bonds are resonating. 

Usually a SCS state is not the exact ground state for 
a given spin Hamiltonian. However if the exact ground 
state is continuously modified into an appropriate SCS 
state, the SCS state describes the essence of the ground 
state. From this viewpoint, the SCS state is useful to 
characterize the phase which the ground state belongs to. 
In some quantum spin chains without side chain, various 
ground-state phases have been successfully described by 
corresponding SCS pictures [19]. Furthermore, a SCS 
state can also quantitatively describe the true ground 
state, if the wave function of each singlet cluster is well 
localized. Among such systems, in the spin system on a 
diamond chain [20], we have the exact tetramer-dimer- 
state solution, which is a kind of SCS state. 



A general SCS picture is defined by a wave function of 
a tensor product form of local singlet states. We call this 
wave function the SCS state and each local singlet state 
a singlet cluster. It is typically written as 

|*)=|t/.(l))® |V(2))---® |^(M)), (29) 

where | '^(p) ) {p = 1, 2, • • • , M) is a singlet cluster and M 
is the total number of singlet clusters in the SCS state. 



C. SCS picture for Gap I phase 

We call the SCS picture for the Gap I phase SCS I. 
The SCS I state is constructed by assuming the following 
requirements: 

(i) The ground state is continuously modified to the 
VBS I state without global rearrangement of the valence 
bond configuration. 
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(ii) The ground state is invariant under the translation 
by a single unit cell. 

(iii) The ground state contains a substantial amount 
of component with J2-valence-bonds in region I-B. 

Assumption (i) is necessary, since the Gap I phase in- 
volves the limiting case of j ^ and k < kc where VBS I 
picture holds, and there is no phase transition in Gap I 
phase. As for assumption (ii), we confirmed that there is 
no indication of the translational symmetry breaking in 
numerical ground states. Assumption (iii) is required to 
explain region I-B. 

Under the above requirements, we take SCS I wave 
function in the tensor product form: 

I \I/i(c) ) =1 Vi(l; c) )® I ^i(2; c) ) • • • ® I MN; c) ) (30) 

as depicted schematically in Fig. 8(a). Here \ipi{p;c) ) is 
a singlet cluster of four | -spins in the pth. unit cell; this 
is denoted by a loop filled in grey. Each singlet cluster is 
a linear combination of two valence bond states written 
as 

|V^i(c))=|ai)+c|A), (31) 

where | ai ) and | f3i ) are the valence bond states defined 
in the right hand side of Fig. 8(b). We have abbreviated 
index p of the unit cell for simplicity. The valence bond 
state I Q^i ) is the same as that shown within a loop in 
VBS I (Fig. 7(a)). The valence bond state \Pi ) contains 
a J2- valence-bond. The coeflicient c should be in the 
limit of j 0, and may be small for region I- A (Ji > J2). 
But it should have a substantial amplitude in region I-B 
(Ji < J2). For finite c, the two valence bond states lo- 
cally resonate in a singlet cluster to contribute to energy 
gain; this effect is examined in subsection VIE, where 
the phase boundary between Gap I and Gap II phases is 
discussed energetically. 

We examined which valence bonds really contribute 
to the ground state by numerically calculating the 
short range correlation functions. For a typical 
case (r = 0.3 and k = 0.2), results are shown in 
Fig. 9. As j increases, |(Si(j9) • S2(p)) | decreases and 
|(Si(p+ 1) • S2(p))| increases. This means that the con- 
tribution from the J2-valence-bonds becomes large in 
comparison with that from the Ji-valence-bonds. It is 
also known that |(Si(p) • T{p -h 1))| takes the maximum 
and |(Si(p) • T(p))| takes the minimum around region I- 
B. Hence a -valence-bond is reduced in region I-B. 
Instead, a valence bond between T{p-\-l) and S]^"'^'^(j9) (or 

S^i\p)) develops, although there is no exchange interac- 
tion between T{p -\- 1) and Si{p). All these results are 
consistent with the SCS I picture. 

D. SCS picture for Gap II phase 

We can the SCS picture for the Gap II phase SCS II. 
Similarly to the case of SCS I, the SCS II state is con- 
structed by assuming the following requirements: 




FIG. 9: Short range correlation functions for r = 0.3 and k 
= 0.2 calculated numerically for N =12, 18, and 24 {N: the 
number of spins). 

(i) The ground state is continuously modified to the 
VBS II state without global rearrangement of the valence 
bond configuration. 

(ii) The ground state is invariant under the translation 
by a single unit cell. 

(iii) The ground state contains a substantial amount 
of component with Ji -valence-bonds in region II-B. 

Under the above requirements, we take SCS II wave 
function in the tensor product form: 

l^ii(c)) =|^ii(l;c))0 |7/;ii(2;c))...0 \iJu{N;c)) 

(32) 

as depicted schematically in Fig. 10(a). Here |'0ii(p) ) is 
a singlet cluster of four |-spins in the pth unit cell; this 
is denoted by a loop filled in grey. Each singlet cluster is 
a linear combination of two valence bond states written 
as 

IV^ii(c)) =|aii)+c|Ai), (33) 

where \ au) and \ f3u) are the two valence bond states 
defined in the right hand side of Fig. 10(b). We have 
abbreviated index p of the unit cell for simplicity. The 
valence bond state \au) is the same as that shown within 
a loop in VBS II (Fig. 7(b)). The valence bond state | l3u ) 
contains a Ji- valence- bond. The coefficient c should be 
in the limit of j ^ 00, and may be small for region II- A 
(Ji < J2). But it should have substantial amplitude in 
region II-B (Ji > J2). For finite c, the two valence bond 
states locally resonate in a singlet cluster to contribute 
to energy gain; this effect is examined in subsection VI E. 

The wave function (33) of the singlet cluster is also 
represented as 

l^ii(c)) = (1 + c) |aii)+c|7ii) (34) 
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(a) 




SCS II '<i^^2 
(Gap II) 

J-\ J2 Si S2 




IPII> 



l«II> 



lYii> 



FIG. 10: (a) SCS II, the SCS picture for Gap II phase (Si = 1, 
& = |, and T = I). SCS II is a tensor product form of singlet 
clusters, (b) A singlet cluster |^ii(c) ) (the left hand side) in 
SCS II. It is represented as a linear combination of two valence 
bond states, | an ) and | Pu ) (the right hand side), c is the 
coefficient of the linear combination, (c) Identity among local 
valence bond states. | pu ) (the left hand side) is exactly the 
sum of I an ) and 1 711 ) (the right hand side) . 



region II-B contains Ji -valence-bonds only in part in the 
linear combination (Eq. (33)). Then, since the energy 
gain owing to Ji -valence-bonds in region I- A is always 
larger than that in region II-B, one might expect that 
region I- A would extend to the whole area of j < 1 . The 
reason why region II-B actually exists is attributed to 
the energy gain by local resonation between the two va- 
lence bond states within each singlet cluster (Eq. (33)). 
Since the resonation is enhanced by the frustration due 
to iir2-interactions, the area of region I-B increases with 
increasing r (= K2/IK1) as really seen in Fig. 2 and 
Fig. 4. If r is fixed, the increase of k (= Ki/2Ji) dimin- 
ishes the effect of Ji -interactions, and makes the effect 
of resonation in SCS II advantageous. This is the reason 
why SCS II appears for relatively large k in the presence 
of frustration r. 



VII. VARIATIONAL CALCULATION 

In the preceding section, we explained that the SCS 
picture represents the essence of each ground-state phase. 
However it does not guarantee that each SCS wave func- 
tion is quantitatively satisfactory. In this section, we 
perform variational calculation using the SCS I and the 
SCS II wave functions to examine the quantitative cor- 
rectness of the wave functions. 



Variational calculation for SCS I 



using identity | Ai ) = | o^ii ) + I 7ii ) 5 which is shown in 
Fig. 10(c). I 7ii ) includes the valence bond on the K2- 
interaction and contributes to energy gain in region II-B 
{K2 > Kt). 



E. Boundary between Gap I and Gap II phases 

First, we consider the j > 1 part of the phase bound- 
ary; it separates regions I-B and II- A as seen in Fig. 6. 
The SCS II wave function in region II-A is close to VBS II 
with J2-valence-bonds, while the SCS I wave function in 
region I-B contains J2 -valence- bonds only in part in the 
linear combination (Eq. (31)). Then, since the energy 
gain owing to J2-valence-bonds in region II-A is always 
larger than that in region I-B, one might expect that re- 
gion II-A would extend to the whole area of j > 1. The 
reason why region I-B actually exists is attributed to the 
energy gain by local resonation between the two valence 
bond states within each singlet cluster (Eq. (31)). When 
k or j increases, the SCS I becomes less favorable and 
VBS II becomes of advantage because of strong Ki- or 
J2-valence-bonds. 

Second, we consider the j < 1 part of the phase bound- 
ary; it separates regions II-B and I- A as seen in Fig. 6. 
The SCS I wave function in region I-A is close to VBS I 
with Ji -valence-bonds, while the SCS II wave function in 



The variational wave function for SCS I is Eq. (30) 
with variational parameter c, which is the coefficient of 
the linear combination. Then the energy per unit cell in 
energy unit Ji is written as 



ei(c) 



1 (^i(c)|i7|^i(c): 

NJi (^i(c)l ^i(c)) 
mc)\h, |V;i(c)) 

{Mc)\Mc)) ' 



(35) 



where H is the total Hamiltonian (1), and h\ is the re- 
duced Hamiltonian for a singlet cluster in SCS I given 
as 



hi = S^') . S2 + iS2 • S"^' + kS\^' ■ T. 



(2) 



j(2) 



(36) 



After straightforward calculation, we have the follow- 
ing formula: 



ei(c) 



3 1 + k - {1 + k + j)c + jc^ 



(37) 



Minimizing ei(c) with respect to c, we have the optimal 
value cn of the coefficient of the linear combination: 



l + k-j-^(l + kf-il + k)j+p 
CI = ^-p^ . (38) 



10 





FIG. 11: Coefficient ci of the linear combination in the varia- 
tional wave function for SCS I {Si = 1, 5*2 = |, and T = |). 
The optimal coefficient c does not depend on r in the wave 
function. The lines are in the order of ascending k from the 
top to the bottom. 



FIG. 12: Coefficient cn of the linear combination in the vari- 
ational wave function for SCS I in the case of r = 0.2 (Si = 1, 
5*2 = |, and T = |). The lines are in the order of ascending 
k from the top to the bottom. 



The result is shown in Fig. 11, where the left axis repre- 
sents — cj. It is independent of r, since a iir2-interaction is 
between different singlet clusters. In the limit of j ^ 0, 
we have ci = and the wave function reduces to VBS I. 
The value |ci| increases with increasing j, and for j ^ 1 
the term of | /^i ) including J2-valence-bonds becomes 
dominant in | ipi (c) ) . This behavior is consistent with 
the argument about SCS I in the preceding section. 



B. Variational calculation for SCS II 

The variational wave function for SCS II is Eq. (32) 
with variational parameter c, which is the coefficient of 
the linear combination. Then the energy per unit cell in 
energy unit Ji is written as 



eii(c) 



1 (^„(c)|i^|^n(c)) 
NJi (^ii(c)l^ii(c)) 
(^11 (c) I hu l^ii(c)) 
(V^ii(c)l V^ii(c)) 



(39) 



where H is the total Hamiltonian (1), and hu is the re- 
duced Hamiltonian for a singlet cluster in SCS II given 
as 



kT . Sf ^ + S^"^ . S2 + jS2 • S^'^ + 2krT • S2. 



2(2) 



(1) 



(40) 



After straightforward calculation, we have the follow- 



ing formula: 



eii(c) = 



3 j + A: - (1 + A; + j - 2kr)c + 

"4 1 - C + C2 



(41) 



Minimizing eii(c) with respect to c, we have the optimal 
value cii of the coefficient of the linear combination: 



cii 



Va2 + a6 + 6^ 
a-\-b 



(42) 



0.2 



with a = j -\-k — l and b = 1 — 2kr. The result for r 
is shown in Fig. 12, where the left axis represents — cn. 
In the limit of j 00, we have cn = and the wave 
function reduces to VBS II. 

The result depends on r, since a i^2-iiiter action is in- 
volved in each singlet cluster. Except for r = 1/2A:, |cii| 
increases with decreasing j, meaning that | Pu) includ- 
ing Ji -valence-bonds becomes dominant in Eq. (33). This 
feature is consistent with the SCS II picture in the pre- 
ceding section. In the special case of 2kr (= K2/J1) = 
1, we have cn = for any values of j, since 6 = in 
Eq. (42). It is explained in the following expression for 
l^ii(c)>: 



|t/'ii(c)) = (l + c)|/3ii)-|7ii; 



(43) 



which is derived from Eq. (33) by identity | Ai ) = | an ) 
+ I 7ii) in Fig. 10(c). For Ji = K2, \ Ai ) including a 
Ji -valence-bond and | 711 ) including a i^2 -valence-bond 
are equally weighted in Eq. (43), as are expected. 
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FIG. 13: Phase boundaries by the variational calculation for 
the SCS wave functions {Si = 1,82 = |, and T = |). 



C. Phase diagram by variational calculation 

In the above subsections, we have obtained variational 
energy ei(ci) of SCS I for Gap I phase and variational 
energy eii(cii) of SCS II for Gap II phase. The phase 
boundary is determined by the equation 



ei(ci) = €ii(cii). 



(44) 



The results for various r are shown in Fig. 13. These 
phase boundaries qualitatively agree with the numerical 
ones in Fig. 4, and support the correctness of the SCS 
pictures. 

Region I-B in the variational phase diagram is wider in 
/c-direction and narrower in j-direction than the accurate 
one (Fig. 4). Further, region II-B is narrower than the 
accurate one. Remembering that regions I-B and II-B 
exist because of local resonation of valence bond states, 
inclusion of longer-range valence bonds resonating with 
each other is possibly effective to improve the variational 
wave functions in regions I-B and II-B. The ground states 
in regions I- A and II- A are close to VBS I and VBS II, re- 
spectively, and the effect of resonation is relatively small. 
Hence, the inclusion of longer longer-range valence bonds 
does not seem to be very effective to improve the varia- 
tional wave functions. 



VIII. SUMMARY AND DISCUSSION 

In this paper, we investigated the nonmagnetic ground 
states of a mixed spin chain with side chains with weak 
frustration. So far, these kind of models have not been 
investigated in depth despite their possible rich physics. 



We have chosen the present model (1) (Fig. 1) as a simple 
and nontrivial one, which will be a good starting point. 
We examined the system in various approaches: a NLSM 
method, a numerical diagonalization method, an inspec- 
tion of limiting cases, an physical interpretation based on 
SCS pictures, and a variational calculation for SCS wave 
functions. 

NLSM methods have been developed for simple spin 
chains without side chain. In the present work, we for- 
mulated a NLSM method for the typical spin chain with 
side chains. The NLSM method analytically provides a 
ground-state phase diagram in the k-j parameter space 
for various values of 6*2, and T. In the special case of 
Si =1, 5^2 = ^, and T = ^, the phase diagram contains 
two quantum disordered phases. Gap I and Gap II, in 
each of which the system has a spin-gap. 

We also examined the case of S^i =1, 5^2 = | and T 
= I by the numerical diagonalization for finite chains. 
Using the method of twisted boundary condition, we have 
determined phase boundaries. When frustration is not 
strong (r < 1), there are two spin-gap phases in the k- 
j parameter space. The numerical results confirms the 
qualitative correctness of the present NLSM method. 

The limiting cases of j ^ and j ^ 00 are precisely 
and analytically treated. For j ^ 1, the Hamiltonian 
(1) describes an array of weakly coupled 3-spin units. As 
k increases from 0, the ground state of each 3-spin unit 
changes from singlet to triplet. Accordingly the whole 
spin chain undergoes a phase transition from the VBS I 
state to the Haldane state. For j ^ 1, the system is de- 
scribed by an effective Hamiltonian where no phase tran- 
sition occurs with changing k as long as r is small. Con- 
sidering the continuity to the large k limit, the ground 
state is described as a state similar to the VBS II state. 

There are regimes where no VBS picture explains the 
ground state for S^i =1, 5*2 = | and T = ^. To explain 
the whole phase diagram, we proposed two SCS pictures; 
SCS I and SCS II for Gap I and Gap II phases, respec- 
tively. Each SCS is a wave function of a tensor prod- 
uct form of singlet clusters. A singlet cluster in both 
the SCS's is a local linear combination of two valence- 
bond states of two different patterns. The resonation 
contributes to the energy gain of the system, and the 
whole phases are consistently explained. 

To quantify the SCS pictures, we performed variational 
calculations with the wave functions representing SCS I 
and SCS II. The phase boundary between Gap I and 
Gap II phases are determined by equating the energy of 
the minimized wave function for SCS I to that for SCS II. 
The resultant phase diagram approximately reproduces 
the phase diagram by the numerical diagonalization. 

Thus we have obtained three phase diagrams: Fig. 2 
by the NLSM method. Fig. 4 by the numerical diagonal- 
ization, and Fig. 13 by the variational calculation. The 
phase diagram of Fig. 4 is accurate, since the extrapo- 
lation by finite size systems is reliable. The other phase 
diagrams qualitatively agree with the accurate one, and 
both the methods are shown to be useful. 
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So far, we have examined the spin chain with side 
chains when the parameter r measuring frustration is not 
large (r < 1). For larger r, wider variety of phases are 
expected. For example, we can extend the analysis in 
the limit of j ^ to the region of r > 1, i. e. kc < 1. 
For kc < k < 1, effective spins of 3-spin units with spin 
magnitude 1 ferromagnetically interacts with the near- 
est neighbors and form a ferromagnetic ground state. In 
terms of the original spins, the ground state is ferrimag- 
netic. Further, by the numerical diagonalization, we have 
found various ferrimagnetic phases with different magne- 
tization in the strongly frustrated regime. The study of 
these ferrimagnetic phases is in progress and will be re- 
ported in a separate paper. 
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